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Abstract. - We study the role of different terms in the A'^-body potential of glass forming 
systems on the critical dynamics near the glass transition. Using a simplified spin model with 
quenched disorder, where the different terms of the real A-body potential are mapped into 
multi-spin interactions, we identified three possible scenarios. For each scenario we introduce a 
"minimal" model representative of the critical glassy dynamics near, both above and below, the 
critical transition line. For each "minimal" model we discuss the low temperature equilibrium 
dynamics. 



S , In the last years many efforts have been devoted to study the relaxation dynamics of 

O ' undercooled liquids near the (structural) glass transition. When the temperature of the liquid 

is lowered down to the critical glass temperature relaxation times becomes exceedingly long and 
difFusional degrees of freedom freeze over very long time scales. As a consequence the difference 
between a structural glass and a disordered system with quenched disorder, which may seem 
;J] ' essential, becomes less and less sharp as the transition is approached since the particles in the 

. . ! liquid become trapped in random position (cage effect) and the dynamics of a single degrees of 

freedom resembles the relaxational dynamics in a random quenched potential. The idea that a 
undercooled liquid is a sort of random solid, which dates back to Maxwell [Q , has been recently 
largely used in the study of the glass transition in undercooled liquids. In this scenario is, for 
example, the study of the Instantaneous Normal Mode (INM) 0, where the N-body potential 
is analyzed in terms of normal modes of oscillations about a given instantaneous configuration. 
Using this technique many properties of the N-body potential in models of undercooled liquids 
have been recently traced out |g, ||, §, y, |§1 . 
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In this letter, based on the connection between structural glasses and quenched disorder 
models, we shall analyze from a general point of view the role of different terms in the iV-body 
potential, such local stress and other non-harmonic terms, on the critical dynamics. We are 
interested only on the general properties therefore we shall use a simple spin-glass model where 
the different terms of the original iV-body potential are mapped into multi-spin interactions. 
In this spirit local stresses are represented by a linear term, the harmonic part becomes a 
two-spin interaction term and higher order nonlinear terms becomes p-spin interactions with 
p > 2. The use of spin models to study the (structural) glass transition, which dates back 
to the end of 80 's 0, has the great advantage that one can construct solvable spin models 
displaying the typical critical behavior of structural glasses [|| |[ |0| . The model we consider 
is a spherical spin-glass model with (random) potential 

p>l \ l<Ji<-<!p<Af 

(v) 

where V^ / , are uncorrelated zero mean Gaussian variables of variance 

/ (,) y^j!pI_ (2) 

and (Ti are N continuous variables obeying the spherical constraint ^^ af = N. The parameters 
Jp define the relative strength of the various terms in (|l|) and can be tuned to reproduce some 
of the properties of the A^-body potential of undercooled liquids obtained for instance by 
INM analysis. This model has been discussed to some extent in literature, see e.g. Refs. 
[ [ll[ [l^ , [l^ , piil , as prototype mean-field model for the structural glass transition, and recently 
also to fit experimental data ||lj]. It is known that the high temperature phase is described 
by the (schematic) Mode Coupling Theory (MCT) for structural glasses. However, to our 
knowledge, a systematic analysis of the relation between the leading terms in the potential 
(M and the critical dynamical behavior, similar to what done in MCT |13], was never done. 
Our analysis identifies three possible scenarios for critical dynamics, similar to what found 
in MCT, and allows for the introduction of a "minimal" spin-glass model for each scenario 
which describes the critical dynamics near, both above and below, the transition. We stress 
that while the high temperature phase is similar to what found in MCT, the low temperature 
phase is different. We finally note that our analysis also reveals all the possible scenarios for 
glassy transition that can be obtained with this model. 

The relaxation dynamics is defined as usual by the Langevin equation 

dta, = -Ra,-^^p^+^,it) (3) 

0(Ti 

where ^i{t) is a Gaussian random field (thermal noise) with variance {^i{t)£,i{t')) — 2 6{t — t'), 
T — 1/(3 the temperature and R a Lagrange multiplier to ensure the spherical constraint which 
must be fixed self-consistcntly. If only the linear term p — 1 and one term with p > 2 are 
present in the expansion (1^) the model is equivalent to the spherical p-spin model introduced 
in Ref . |q| whose dynamics has been studied in Ref. 0, Ha . The case where only the quadratic 
term p = 2 is present has been analyzed in Refs. |l7l ]lq|. It is worth noting that there is a 
nontrivial difference between the dynamics for p = 2 and p > 2. Indeed while in the former 
case relaxation is an orientation process of the state vector {ai} toward a (doubly degenerate) 
state p7| , for p > 2, similar to what happens in structural glasses, relaxation occurs in a 
free energy landscape characterized by many (highly degenerate) local minima M, M . In this 
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spirit small non- harmonic terms change qualitatively the dynamical properties, therefore in 
this letter we always assume at least cubic nonlinearities. 

The analysis of (0) simplifies considerably in the thermodynamic limit iV — > oo, where the 
dynamics can be described by a set of self-consistent equations involving a single spin only 
and the averaged correlation and response functions C{t, t') and G{t, t') M, M, 0, pl[ . 

In the high temperature phase at equilibrium C and G are related by the fluctuation 
dissipation theorem (FDT) G{t,t') = G{t - t') ^ -0{t - t')dtC{t - t') and the the mean-field 
dynamical equation simplifies further, for details see e.g. Refs. [g^, 0, |^, ^. The resulting 
mean-field dynamical equation can be written as, 

[dt + r] G{t) + [ ds A{t - s) dsC{s) = 1 - r (4) 

Jq 

where A{t) = A[C(t)] = Ep>2MpC^~'(i), r = R- Ep>2t^P^ Mp = JpP/2T^, and t is 
now a time difference. For large t we can define the Edward-Anderson order parameter as 
limt^oo G(t) — go which, taking the t — > cxd limit of eq. (]^, obeys the equation: 

Ml + Afeo] = 7j — ^^^ (5) 

The parameter r has been eliminated using the the spherical constraint which now reads 
C(0) = 1. This is the "replica symmetric" solution for this model. Stability analysis reveals 
that this solution is stable iff 

dMqol < _J__ (g) 

dqo (1 - qoy 

the equality being satisfied along the transition line. Since the linear term in (|l|) acts as an 
external field, it can be shown that if /ii = the only stable solution is go = while for /ii > 
we have < q'o < 1 [§| • The order parameter go is the time-persistent part of the correlation 
induced by the variance of the local stress, therefore go 7^ is not associated to a glassy phase. 
We note that in structural glasses has been found [Q that local stresses develop a zero mean 
Gaussian distribution approaching the glass transition, therefore even if their role is irrelevant 
well inside the liquid phase, they can be relevant for the dynamical transition. To study the 
relaxation near the transition we introduce the (rescaled) connected correlation function (j>{t) 
by writing C{t) = go + (1 — qa)4>{t) where, from eq. (Q) and (||), (f>{t) obeys the equation 

dt(j){t) + (f){t) + / ds AA{t - s) ds(l>{s) = (7) 

1-90 Jo 

with 0(0) = 1 and 

oo 

AA(t) = A[G{t)] - A[qo] = 5] A,/ it) (8) 

fe=i 

This equation has the same structure of the schematic MCT equation for glasses considered 
by Gotze pSl. The factor (1 — go)^^, absent in the Gotze equation, just changes the form of 
the long time solution near the type B transition (see below). The Gotze equation is recovered 
for go = 0. 
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Fig. 1. - Schematic time behavior of the correlation hmction C{t) near the type A transition 
(a) and the type B transition (6)). 



In deriving eq. (m) we made use of FDT so that this equation is appropriate only above 
the dynamical transition. In the low temperature phase (spin-glass or glass phase) the FDT 
must be modified M, n& and, moreover, non-equilibrium p6[ and equilibrium [Bl dynamics are 
separated by infinite time scales and described using different approaches. 

The transition from the high to the low temperature phase is of type A or B depending on 
the instability of eq. (^ . The type A transition occurs along the instability line of solution 
(ph [equal sign in eq. (0)] and is the analogous of the De Almeida and Thouless line in spin 
glasses. Type B transition appears with a sudden appearance of another long time persistent 
solution: limt_^oo '/'(i) = / 7^ where, from eq. (0), / is solution of the bifurcation equation: 
AA[/] = / [(1 - (7o)(l - f)]~^- For go = this reduces to that found by Gotze (23|. Defining 
/ — {qi —qoVl — qo)^^ we recover the usual equation for p-spin-like models at the discontinuous 
transition |y, ^. The relaxation dynamics near the two types of transitions is quite different, 
as shown schematically in Fig. 1 where relaxation parameters are defined. 

Which kind of transition - type A or B - takes place depends on the parameters fip. Near 
the transition only the first two non-zero terms in the sum (M are relevant, all others giving 
sub-leading corrections which can be neglected at the transition. The six non-trivial cases are 
reported in table |. All other cases can be qualitatively mapped to one of these. With /Xp>3 
we mean a generic non-harmonic term of order p > 3. The six cases can be grouped into three 
different classes depending on the two relevant A^. Following a notation introduced by Gotze, 
we call them 1 — 2, 1 — 3 and 2 — 3, respectively. For each class we can define a "minimal" 
spin-glass model obtained by retaing only two terms in the sum (n|) . If we denote these modes 
with "x + y"-SG model, where p — x,y are the two retained terms, the simplest choice for 
spin-glasses "minimal" model is: 2 -|- 3-SG model for class 1 — 2, 2 + 4-SG model for class 1 — 3 
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Table I. - = missing, 1 = present , - — irrelevant 



and 3 + 4-SG model for class 2 — 3 [table ||. With this choice the relation between /I's and 
A's is very simple: A^ = /ifc+i- One could use other choices, e.g. the 1 + 3-SG model for the 
class 1 — 2. This would lead to a different relation between fi's and A's, but will not change 
the results if expressed in terms of A's. 

The "minimal" models can be analyzed in full details also far from the transition and hence 
can be used to gather more informations on the low temperature phase in the different classes. 
For this reason in what follow we shall restrict to the "minimal" models and consider the 
equilibrium dynamics in the low temperature phase. In is clear that the results we shall obtain 
are also valid for the full model as far as the extra terms in potential Im) can be neglected, 
i.e., only near the transition. The equilibrium dynamics in the low temperature phase can be 
studied using the technique of Ref. pOl H| which assume equilibration on very long times and 
a modified form of FDT for the very slow processes. This method is known to reproduce the 
low temperature phase statics. The calculation is lengthly and will not be reported. Here 
we quote the main results. The non-equilibrium aging phenomena can be studied using the 
technique of Ref. [|l6|, but these are beyond the purpose of this paper. 

In Fig. 2 (a) we report the phase diagram in the Ai — A2 plane of the "minimal" 2 + 3-SG 
model. This is equal to the phase diagram of the spherical p-spin model in a field 0, M when 
expressed in terms of Ai_2- For small A2 a transition of type A separates a paramagnetic 
from a glassy phase described by "1-step" replica symmetry broken solution. By increasing A2 
along the critical line Ai = 1 we eventually reach a tricritical point where the type B solution 
appears and the transition changes from A to B. Since / is zero at this point there is no jump 
in the bifurcation parameter. As far as the transition lines are concerned, this phase diagram 
is similar to what obtained in MCT for structural glasses |23] , however the properties of the 
low temperature phase are different M, 0]. 

The phase diagram of the "minimal" 2+4-SG model, shown in Fig. 2 (b), is more reach. The 
high temperature phase is still separated from the low temperature phase by a type A transition 
for small A3 and by a type B transition for higher A3. However the low temperature phase 
consists now of two different glassy phases: one described by "1-step" and one by "00-steps" 
replica symmetry broken solution. The transition from 1-step to 00-steps RSB appears as an 
instability of the 1-step RSB phase - type A - and in this respect is similar to what found 
for the Ising p-spin model p3. The point where the type A and type B transitions meet is 
not a critical point since the two lines just cross, and / jumps discontinuously. The possibility 
of two different glassy phases in model (P) was discussed in Ref. |13] from considerations on 
the statics of the model, but the transition line between the phases was not determined. We 
stress that the observed low temperature phase is completely different from what predicted 
by Gotze and Sjogren |n5[ using the MCT equation M). Indeed they find two "1-step" RSB 
phases separated by a type B transition line ending at a critical point. This, however, does 
not contradicts our results since as previously discussed in the low temperature phase eq. (Q) 
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Fig. 2. - The phase diagram of the "minimal" 2 + 3-SG model (a), and of the "minimal" 
2 + 4-SG model (6). The exponents of C{t) [see Fig. 1] and the value of / are shown at the 
edges of the transition lines. 



is not appropriate for our model. 

The "minimal" 3 + 4-SG model is the simpler one since the absence of quadratic term in 
the potential makes the continuous transition impossible K^. A discontinuous transition of 
type B separates the high temperature phase form a low temperature phase of "1-step" RSB 
type, see Fig. 3. 

In conclusion in this Letter we have shown that simple multi-spin interactions spin-glass 
spherical models, which obeys MCT equations above the dynamical transitions, lead to three 
different scenarios for critical dynamics according to the leading terms in the potentials. The 
type of dynamical transitions can be continuous as well as discontinuous. From a general point 
of view a continuous transition may occur only if harmonic and/or linear terms are important. 
Within the assumption that different terms of the real iV-body potential of glass-forming 
systems can be mapped into the multi-spin interaction terms in (||), using the informations 
from INM studies the type of scenario for the critical dynamics can be predicted. 

Schematic models of MCT for structural fragile glasses and undercooled liquids belong to 
the 1 — 2 class of universality. This is in agreement with our conclusions since in these systems 
local stresses and cubic nonlinearities are known to be important B. Based on the same 
assumptions the continuous transition observed in the rotational dynamics of linear molecules 
[ p5| p^ could be related to the relevance of linear and/or harmonic terms in the rotational 
degrees of freedom potential. To test our conjectures it would be also of interest to find physical 
systems belonging to 1 — 3 or 2 — 3 dynamical universality class. 

To get more insight on the low temperature phase for each scenario we have introduced a 
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Fig. 3. - The phase diagram of the "minimal" 3 + 4-SG modeL The exponents of C{t) [see Fig 
1] and the value of / are shown at the edges of the transition lines. 



"minimal" SG model and discussed the equilibrium low temperature dynamics. In particular 
for the 1 — 3 class has been calculated the critical transition line between the two glassy 
phases. The results found for the "minimal" models do apply to the full model near the 
transition lines, both above and below the transition. Recent results p7| have shown that 
finite-size mean-field p-spin-like models, where activated processes are allowed, exhibit strong 
similarities with structural fragile glasses. Therefore the "minimal" models can be highly 
valuable to study the glass transition in the different scenarios also beyond MCT. Work in this 
direction is in progress. 

We finally stress that there are examples of glass forming systems where the MCT scenario 
is different form those proposed in the present Letter. For example a mixture of sticky 
hard-sphere exhibits a two different glass phase separated by a type B transition [ P8| . In 
this case the short wave-length dependence of vertices in the MCT plays a crucial role and 
therefore a simple model where inhomogeneous spatial fluctuations are neglected, as the one 
considered in this Letter, cannot be used. 
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